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ABSTRACT 

Observations of an optical source coincident with gravitational wave emission detected 
from a binary neutron star coalescence will improve the confidence of detection, pro- 
vide host galaxy localisation, and test models for the progenitors of short gamma 
ray bursts. We employ optical observations of three short gamma ray bursts, 050724, 
050709, 051221, to estimate the detection rate of a coordinated optical and gravi- 
tational wave search of neutron star mergers. Model i?-band optical afterglow light 
curves of these bursts that include a jet-break are extrapolated for these sources at 
the sensitivity horizon of an Advanced LIGO/ Virgo network. Using optical sensitiv- 
ity limits of three telescopes, namely TAROT (m=18), Zadko (m=21) and an (8-10) 
meter class telescope (m=26), we approximate detection rates and cadence times for 
imaging. We find a median coincident detection rate of 4 yr _1 for the three bursts. 
GRB 050724 like bursts, with wide opening jet angles, offer the most optimistic rate 
of 13 coincident detections yr _1 , and would be detectable by Zadko up to five days 
after the trigger. Late time imaging to m — 26 could detect off-axis afterglows for 
GRB 051221 like bursts several months after the trigger. For a broad distribution of 
beaming angles, the optimal strategy for identifying the optical emissions triggered 
by gravitational wave detectors is rapid response searches with robotic telescopes fol- 
lowed by deeper imaging at later times if an afterglow is not detected within several 
days of the trigger. 

Key words: stars - gamma-ray burst: individual - gravitational waves - techniques: 
miscellaneous - stars: neutron 



1 INTRODUCTION 

A multi-messenger approach to gravitational wave detection 
is one of the most prioritized goals for second generation 
ground-based gra vitational wave (G W) detectors, such as 
Advanced LIGO (jHarrv et al.1 12010T ) and Advanced Virgo 
jAcernese et alj[2009f ) . It allows GW candidates that are too 



weak to claim detection based on GW data alone to be asso- 
ciated wit h an optical signal that cou ld provide strong con- 
firmation. iKochanek Sz Piran I (| 19931 ) showed that in prin- 
ciple a joint electromagnetic-gravitational wave (EM-GW) 
search reduces the GW amplitude detection threshold by 
a factor of about 1.5. This has the effect of extending the 
sensitivity horizon distance of GW detectors for a binary 
neutron star merger (NS-NS), so that the number of poten- 
tially detectable GW sources increases by a factor of ~ 3.4. 
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The benefits of joint EM-GW searches are significant on 
many fronts. For example, direct measurement of GWs from 
a NS-NS late-inspiral and merger provide a means to deter- 
mine luminosity distance; by combining GW derived lumi- 
nosity distance measurements with EM redshifts, one can 
constrain key cosmological parameters. Joint observations 
offer unprecedented insight into the complex astrophysics 
that are key to the EM emissions, offering the most com- 
plete picture from the strong to weak gravity regime. 

One expected EM counterpart of a NS-NS merger is 
a short gamma ray burst (SGRB), where 'short' is defined 
as T90 < 2fl The popular model for SGRBs is a compact 
object merger tri ggering an explosion causing a burst of 
collimated 7-ra ys (jEichler et al 1 1l989l : iNaravan et all 1 19921 : 
iLee et al"1l2005l ) powered by accretion onto the newly formed 
compact object. SGRBs are believed to be produced by dis- 
sipation of kinetic energy of ultra-relativistic outflow from 
the central engine with a Lorentz factor of T ~100-1000. 
The outflow is eventually decelerated by interaction with 
interstellar matter to produce a fading x-ray and optical af- 
terglow. After r decreases to T ~ where 6j is the jet 
opening half angle, the radiation beam is wider than the 
outflow, so the afterglow becomes observable from angles 
greater than Qj. 

Although NS-NS mergers are the favoured progenitor 
for SGR Bs, other scenarios cannot b e excluded, i.e. NS-BH 
mergers (IMcWilliams fe Levinll201ll ) or magnetar outbursts 
(Nakar 2007). The main evidence is based on the associ- 
ation of some SGRBs with an older stellar population, as 
compared to 'long' GRBs, which are associated with mas- 
sive stellar collapse. Evidence for the origin of SGRBs in the 
final merger stage comes both from the host galaxy types 
(e.g. Zheng & Ramirez-Ruiz 2007) and the measured off- 
sets of GRBs from their host galaxies (e.g. Belczynski et al. 
2006) . Kicks imparted to NSs at birth will produce velocities 
of several hundred km s~ , implying that binary inspiraling 
systems may occur far from their site of origin. Fong, Berger 
& Fox (2009) using Hubble Space Telescope observations to 
measure SGRB galaxy offsets, find the offset distribution 
compares favorably with the predicted distribution for NS- 
NS binaries. However, they do not rule out at least a partial 
contribution from other progenitors systems. 



1.1 GW triggered SGRB afterglow search 

A number of laser interferometric GW detectors have now 
reached their design sensitivities and have been operating 
as a global array, coordinating with electromagnetic ob- 
servations through triggered follow-ups. These include the 
LIGC0 detectors based at Hanford and Livingston in USA, 
the Virgrjf] detector in Italy and the GEO 6OC0 detector in 
Germany. The LIGO and Virgo detectors are undergoing a 
series of upgrades towards Advanced configurations that will 
produce an order of magnitude improvement in sensitivity. 



Advanced LIGC0 and Advanced Virgr]f| are expected to be 
operational by 2015. 

Because the coalescing binary NSs are expected to radi- 
ate GWs in the sensitivity band of Advanced LIGO/ Virgo, 
coincident GW-EM observations of SGRBs will determine 
if the engine is a NS-NS or NS-BH binary merger. Fur- 
thermore, the rates of EM and coincident GW detections 
could constrain the distribution of jet collima tion angles of 
SGR Bs, crucial for understanding energetics (|Berger et al I 
l2007f ). This is possible because the binary inclination an- 
gle to the line of sight is a GW observablfl A direct con- 
sequence of collimation is that the rate of (both long and 
short) GRB afterglows should be higher than those ob- 
served as prompt bursts. SGRBs afterglows observed 'off- 
axis' without a prompt counterpart have been termed 'or- 
phan afterglows'. There has not been a definitive discovery 
of an orphan afterglow, despite both dedicated searches to 
m = 23 and using archived data (|Totani fe Panaitescu|[2002l ; 
iRau. Greiner fe Schwarzl 120061 ). The non-detection could 
partly be attributed to the small flux of off-axis afterglows, 
compared to on-axis ones that can be observed as early as 
seconds after the prompt burst and in some cas es while the 
high- energy prompt emission is still occurring (|Klotz et al.l 
120091 ). 

The prospects might seem bleak for detecting off-axis 
afterglows. This is the case for an all-sky search, but a GW 
observation of a NS-NS inspiral collapses the search area 
from all-sky to the error ellipse (angular sky resolution) of 
the GW detectors. For such events, the angular resolution 
of a detector network depends on GW source strength, ori- 
entation of the binary axis, and on the ge ometrical configu- 
ration of the network |Wen fe Chen|[201ol ). Error ellipses of 
the order of tens to a few square degrees can be achieved 
through triangulation of time differences i n signal arrival 
times at various det ectors in a network dGiirsel fc Tintol 
ll989l ; lFairhurstll2009l ). 

The sensitivity distance, D s , for sources uniformly dis- 
tributed in orientation and sky location is approximated by 
-Dh/2.26, were Dh is the horizon distance in Mpc at which 
an optimally-orientated, overhead source can be det ected 
with a signal-to-noise ratio of 8 l|Abadie et al.ll2010l ). We 
note that although Dh is based on the sensitivity of sin- 
gle detector, non-optimal effects such as non-Gaussian, non- 
stationary detector noise allow the approxima tion to be as- 
sume d for an Advanced LIGO /Virgo network (|Abadie et alj 
l2010t ). Significantly, the value D s of Advanced LIGO/ Virgo 
for a NS-NS coalescence is smaller than the average distance 
to the observed SGRB population. This implies that EM 
emissions from NS-NS coalescences triggered by Advanced 
LIGO /Vir go will be brighter t han t he triggered GRB 

emissions. iNuttall fc Suttonl (|2010l ) showed that a galaxy 
ranking procedure could identify the host galaxy 75-95% of 
the time out to 100 Mpc for 5 images taken by narrow field 
and wide field telescopes respectively. Their method depends 
on galaxy survey completeness, so may not be applicable for 



1 The duration in which the cumulative counts increase from 5% 
to 95% above background. 

2 http://www.ligo.caltech.edu/ 

3 http://www.virgo.infn.it/ 

4 www.geo600.uni-hannover.de 



5 www.ligo.caltech.edu/advLIGO/ 
B www.cascina.virgo.infn.it/advirgo/ 

7 In practice, for sources not associated with host galaxies, the 
inclination angle has a strong degeneracy with distance, particu- 
larly for angles less than 45°. 

8 http:/ /swift. gsfc.nasa.gov/docs/swift/swiftsc. html 
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the Advanced GW detector network searches. We discuss the 
issues that will affect the localisation of SGRBs in Section 
II 

To estimate the optical flux of a SGRB as the EM coun- 
terpart o f a NS-NS coalesce nce, we use the plausible esti- 
mates of (jAbadie et alJl2010h for the Advanced LIGO/Virgo 
sensitivity distances and detection rates of NS-NS coales- 
cences. Taking a rate density of NS-NS coalescences ~ 
10" 6 Mpc- 3 yr _1 and D H = 445 Mpc, they find D s » 200 
Mpc and a detection rate Rdet ~ 40 yr _1 . This rate could 
potentially be increased by considering the improved signal 
to noise ratio for a coincident GW and optical search. The 
estimated increase in signal to noise ratio is about 1.5, as- 
suming a narrow coincidence window, but the optical after- 
glows may not be imaged until hours after the GW trigger. 
Nonetheless, the afterglows should be relatively bright at the 
distances we are considering so it is possible that light curves 
could be extracted and extrapolated to earlier times. Hence, 
we assume the sensitivity distance increases by a factor of 
1.5, to 300 Mpc, so that 7? dc t ~ 135 yr" 1 . 

In this Letter we investigate the temporally varying op- 
tical brightness of SGRBs at the Advanced LIGO /Virgo sen- 
sitivity distance to NS-NS coalescences, from the first hour 
of the EM emission to about a hundred days later. We use 
optical data for three SGRBs, GRB 050724, GRB 050709 
and GRB 051221, that show evidence for collimated emis- 
sion and are localised. Assuming that all NS-NS mergers 
produce jets, we use these data combined with reasonable 
estimates of the NS-NS coalescence rate detected by a GW 
detector network to determine the fraction of events that 
could be potentially detected as both on and off axis bursts. 
For definiteness, we use the sensitiv ity limits of two r obotic 



18, see iKlotz et all liooi ). and 



telescopes, TAROT (m 

Zadko (m = 21, see lCoward et al. II2010T ). that have partic- 
ipated in the optical follow-up of LIGO/Virgo GW triggers 
as pa rt of the LOOC UP program during the 6th Science 
Run ijKanner at al.ll2008h . We also include a much deeper 
sensitivity limit of m = 26, representative of an (8 — 10)m 
class telescope such as the VLI0. Finally we discuss the 
main issues that need addressing to interpret and optimize 
the science return from joint optical and GW searches. 

2 OPTICAL AFTERGLOW MODEL AND 
COINCIDENT RATES 

2.1 Observations 

In order to constrain the detection rate, we require locali- 
sation (including redshift), beaming angles and the optical 
flux values for our sample of bursts. We use three SGRBs 
that have estimates for these parameters: namely GRB 
050709, GRB 050724, and GRB 051221A. 

GRB 050709: From comparison of X-ray and optical 
data, a jet break is present in the optical at about 10 
days after the burst (Fox et al. 2005). On the other hand, 
Watson et al. (2006) claimed that the light curves were 
not displaying such a break. We note however that they 
excluded one optical data point within their fit, arguing it 



was coincident with a late X-ray flare. However, the data 
point they excluded was 9.8 days after the burst, compared 
to 16 days for the X-ray flare. We assume the explanation of 
Fox et al. (2005), noting that the detection of the jet-break 
is supported by only one data point. 

GRB 050724: From radio and near infrared data, 
Panaitescu (2006) and Berger et al. (2005), claimed evi- 
dence of a jet break about 1 day after the burst. However, 
the main feature of the jet-break is its achromacity (Rhoads 
1999), and the X-ray data do not feature such a break 
(Gruppe et al. 2006). The X-ray light curve is consistent 
with no jet break up to 22 days after the event, or 8j > 25°. 
GRB 051221A: The detection of the jet break was 
observed in X-ray only (Soderberg et al. 2006). A jet-break 
is clearly visible in the light curve at about 5 days post-burst. 

We note that determination of the jet opening angle 
from the break-time is strongly sensitive on the model pa- 
rameters used. In the previous cases, the model was al- 
ways the same (a forward shock fireball expanding with ki- 
netic energy E in a medium of constant density n), where 
6 oc (n/E) 1 ^. However, n and E are uncertain in all bursts 
(see e.g. Panaitescu 2006 for the values of n in the cases of 
GRB 050709 and GRB 050724). Hence, opening angle val- 
ues reported in these papers are taken as indicative values, 
and we use them as a guide (this is why we do not perform 
a k-correction for cosmological effects, as the correction is 
small compared to the uncertainties of the beaming angles). 

2.2 Afterglow model and rates 

For GRB 050724, GRB 050709 and GRB 051221, 
we define Fiihr using _R-band flux es at 11 hr from 
iNvsewander. Fruchter fc Pe'er I (|2009T) . but scaled to a 
source distance of 300 Mpc (the mean sensitivity distance of 
the GW search). Equation (1) scales the flux using a power 
law index a± until jet-break time t ^ tj t o, when the beamed 
emission starts to expand sideways, and by the luminosity 
distance to the source dh(z): 



D a 



(1) 



The above equation can be used to estimate the _R-band 
flux at times before the jet break. To model the SGRB light 
curve at post jet-break times, we employ a smoothly joined 
broken power law (Beuermann et al. 1999), 



F(t) = Fj 



+ 



-i/t. 



(2) 
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where Fj is the flux at the jet break time tj, a 1 and o<2 are 
the pre-break and post-break light-curve slopes, and n scales 
the sharpness of the break. The beaming angles and break 
times for GRB 050724, GRB 050709, GRB 051221 are (25°, 
22d), (14°, lOd), (7°, 5d) respectively with corresponding 
power law indices of ai and Q2 ( — 1.5, —2), ( — 1.25, —2.83) 
and (-1,-2). 

Table 1 shows the derived parameters with the extrap- 
olated _R-band magnitude at 1 hr post burst at 300 Mpc. 
We point out that the optical data used in the references to 
derive the beaming angle and break times is uncertain, and 
we do not account for optical bumps and flares that can be 
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Table 1. The main observed and derived parameters of GRB 
050724, GRB 050709 and GRB 051221. See SectionOfor caveats 
and uncertainties in Q\ . 



GRB 


log E iao 


z 


R-mag (1 hr)t 


tj (d) 




050724 


50.21 


0.26 


12.7 


< 22 


25 


050709 


49.06 


0.16 


17.2 


10 


14 


051221 


50.95 


0.55 


13.7 


5 


7 



t magnitudes are converted from flux (Jy) to the AB magnitude 
system using hiab = -2.5 log(F) + 8.9 at a source distance of 



300 Mpc 



significant, especially at early times. Nonetheless, it is clear 
from Table 1 that if one of the well localised SGRBs occurred 
within D s , and was on-axis, it would be bright at early times 
and easily detected by modest aperture telescopes. 

To calculate the rate of triggered detection of on-axis af- 
terglows requires accounting for the beaming angle 9j . Equa- 
tion © calculates the number of afterglows seen per yr 
for a certain 9j assuming a LIGO/ Virgo detection rate of 
Rdet ~ 135yr _1 (see section [TTTj) 



-Ron = -Rdct[l - COs(6>j)] . 



(3) 



The rate of off-axis bursts is determined by the maxi- 
mum angle, # max , away from the jet center that an off-axis 
observer could see the afterglow. This angle depends criti- 
cally on the flux limit of the telescope Fy lul , Fj and 9j. For 
off-axis detection, the main constraint is the off-axis emis- 
sion only becomes visible after the j et-break time, which 
may b e days after the prompt burst. iTotani fc Panaitescul 
(2002), assuming the uniform jet model, show that 6 m « can 
be expressed as 



-(3+5) 
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Fun 



-1/(202) 



(4) 



where «2 is the post-break optical decay index and 8 is a nu- 
merical factor ~ 1. 0j can be replaced with (9 max in equation 
Q to estimate the SGRB off-axis detection rate; 



Ro 



= #dct[l - CQS( 



R 



(5) 



We note that if # max ^ 8j only on-axis afterglows will be 
visible at any time, and the rate is determined by 0j. 

3 RESULTS 

Using the above relations and light curve characteristics for 
GRB 050724, GRB 050709 and GRB 051221, we extrapolate 
the light curves beyond the jet-break times to constrain de- 
tection limits, rates and cadence times using the sensitivities 
of TAROT, Zadko and an (8 — 10)m class telescope. Figure 
1 shows the temporal evolution of the three R — band light 
curves using equations (1) and (2) at a source distance of 300 
Mpc, and published values for the decay indices. The three 
curves are quite different, GRB 050724 and GRB 051221, 
are relatively bright at early times, and can be seen from 
days to some tens of days by meter class telescopes. 

Table 2 shows # max , the detection rates R on and Rag, 
and the maximum possible times that an on-axis burst 
would be visible for the three telescopes. We find that GRB 
050724 like bursts are detectable at a relatively high rate. 
This is because the initial beaming angle is large > 25°, so 
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Figure 1. Three model light curves using equation (0 for 
GRB 050724, GRB 050709 and GRB 051221 extrapolated to a 
source distance of 300 Mpc, the horizon limit for the Advanced 
LIGO/Virgo detector network. The beaming angles and break 
times for the model bursts are (25°, 22d), (14°, lOd), (7°, 5d) 
respectively. Power law indices before and after the breaks are 
(-1.5, -2), (-1.25, -2.83) and (-1, -2) respectively. The hori- 
zontal dashed lines from bottom to top are the approximate sen- 
sitivities for an (8-10)m class telescope, Zadko Telescope (lm) 
and TAROT (0.25m) respectively. The maximum on-axis times 
for visibility are shown in Table [2] 



more of the afterglows can be seen on-axis. Conversely, the 
brightest afterglow at the post break time, GRB 051221, is 
the least likely to be detected because of the small 8j ~ 7°. 
It is apparent that both TAROT and Zadko are unlikely 
to detect off-axis afterglows from these SGRBs, but a tele- 
scope capable of deep imaging to m = 26 could detect an 
additional 5 afterglows yr -1 for GRB 051221 like events. 

Table 2 also shows the maximum time, t max , that the 
telescopes could detect the SGRB afterglows. This sets the 
limit on the cadence times for imaging. GRB 051221 has the 
brightest afterglow and is potentially detectable the longest 
time; t max ~ lid for Zadko. Given that the GW error ellipse 
is of order degrees in size, identification of a transient is more 
feasible for this afterglow type. Unfortunately, they occur 
at a rate of 1 yr -1 , and given optical selection effects (see 
section [4] for a discussion) may be missed altogether. GRB 
050724 like events, occurring at an optimistic rate of 13 yr" 1 , 
would be detectable up to 5d by Zadko. This would allow 
time for surveying degree size fields and multiple telescopes 
at different longitudes to perform follow-up imaging. 

Our results in Table 2, based on a broad distribution 
of beaming angles, suggest that the optimal strategy for 
identifying the optical emissions triggered by gravitational 
wave detectors is through initial rapid response searches 
with robotic telescopes, followed by deeper imaging at later 
times if an afterglow is not detected within several days of 
the trigger. 

4 DISCUSSION 

The first attempts for a triggered search of the optical coun- 
terparts of NS-NS coalescences using GW detectors are just 
commencing. There are many uncertainties and issues that 
will need careful consideration for these types of searches. 
Firstly, our results show that the coincident detection rate 
depends critically on the beaming angle distribution. For 
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Table 2. Equations (4) is used to calculate the maximum off- 
axis angle that the three SGRBs could be observed assuming 
limiting magnitudes of TAROT, Zadko and an (8 — 10)m class 
telescope such as VLT. Equations (3) and (5) are used to calculate 
both on and off-axis rates of SGRB afterglows associated with 
GW emission form NS-NS coalescences detected by a 3-detector 
network of GW detectors. The maximum time, t max for detecting 
the optical afterglows for the three telescope sensitivities is also 
shown. The median and most optimistic on-axis detection rates 
are 4 and 13 yr — 1 respectively. 
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TAROT 
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1 





2 
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11 




VLT 


18 
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The maximum off-axis viewing angle, 6> max , is not shown for 

afterglows with max < 8j. This also implies a non-detection 
of an off-axis afterglow. The rates are upper limits and do not 
account for Galactic extinction and crowded fields. 

nearly isotropic optical emission, similar to GRB 050724, 
the coincident rates are very promising and will improve the 
confidence of the GW detection and provide much needed lo- 
calisation. Even non-detections of optical emissions for high 
signal-to-noise ratio NS-NS GW candidates is interesting. 
Non-detections of a statistically significant sample would 
constrain the SGRB beaming angle distribution, or show 
that SGRBs are not linked to NS-NS coalescences. Both 
implications are critical to our understanding of the progen- 
itors of SGRBs. 

To fully use joint optical and GW observations requires 
understanding and accounting for selection effects that have 
historically plagued SGRB optical observations. Typical se- 
lection effects include Galactic extinction and crowded star- 
fields. There is also the problem of host galaxy extinction, 
although this is expected to be a less of a problem given 
that at least some of the afterglows are observed off-set from 
the hosts. From the current optical follow-up attempts of 
Swift triggered bursts, it is clear that a significant fraction 
of SGRBs that are apparently on-axis have not been ob- 
served in the optical at all. This can partly be attributed to 
the much greater distances of the Swift bursts, compared to 
a GW triggered search. We have based our calculations on 
a sensitivity distance of ~ 200 Mpc (non-coincidence rate), 
but note that GW wave sources with smaller angles of in- 
clination will be detected at greater distances. Preliminary 
studies suggest that this bias can increase the median dis- 
tance to ~ 269 Mpc for binaries with inclination angles less 
than 25°. A numerical study is planned to investigate how 
this bias will effect the rate of coincident detections. 

Another important issue for the joint searches is the 
large errors in the GW source localisation, which can extend 
to some tens of degrees. This is a significant problem given 
the first observed GRB optical afterglows required small er- 



ror boxes of arc-minute size. To counter this, the GW trig- 
gered search strategy will use the estimated horizon distance 
of the detector network to reduce the number of potential 
host galaxies, as opposed to a 'blind' error box that extends 
to cosmological distances. This technique was demonstrated 
in the recent LIGO/Virgo S6; the telescope pointing strategy 
used a galaxy weighting taking into account the mass and 
the distance of catalogued galaxies. A detailed description 
and analysis of the S6 coincidence searches will be published 
as a LIGO Scientific Collaboration paper. 

Another problem that manifests with large coincidence 
error boxes is the increasing chance of detecting false coinci- 
dent optical transients. False coincident sources may include 
supernovae, flare stars, variable active galactic nuclei and 
even Earth orbiting space debris. Fortunately, some of these 
sources can be excluded in the analysis because of the sen- 
sitivity distance of GW searches and the expectation that 
the strongest GW sources will be associated with catalogued 
host galaxies. Another possibility for optimizing the coinci- 
dence search is to data mine archived images from very wide 
field optical surveys, such as SkvMavver\ ll \ and into the fu- 
ture the planned Large Synoptic Survey TelescoveV\ 

Given the above practical difficulties, an observed as- 
sociation of an optical transient with a NS-NS coalescence 
triggered by Advanced LIGO/Virgo is challenging. But, the 
science pay-off for such a discovery is enormous and pro- 
vides motivation to address the issues discussed here in more 
detail. Now is the time to determine the optimal strate- 
gies for optical follow-up in readiness for the more sensitive 
GW searches in the following years. To accomplish this will 
require a more comprehensive understanding of optical se- 
lection effects, the false alarm rate expected from SGRBs 
within the error ellipses of GW networks, and techniques to 
improve the localisation of the host galaxy. 
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